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1. Introduction 

Since the original work of Moore and Shannon p], network reliability has been a 
field devoted to the calculation of the connection probability between different sites 
of a network constituted by edges (links, bonds) and nodes (vertices, sites), each of 
them having a probability of operating correctly (the reliability). This field, although 
mainly developed in an applied background |2|, is strongly related to graph theory 
0, H], combinatorics and algebraic structures [5j E], percolation theory [3, [8], as well 
as numerous lattice models in statistical physics 0, [TOj [TTJ [12]. For instance, the all- 
terminal reliability ReLi, i.e., the probability that all nodes are connected, is derived 
from the Tutte polynomial, an invariant of the associated graph, when all edges have 
the same reliability p (0 < p < 1). This polynomial appears in the partition function for 
various Potts models, and has been calculated for several families of graphs j9j [TO, EH]; 
the location of its complex zeros has also been studied [TQl HH EEH]- The two-terminal 
reliability Rek^s — > t), the probability that a source s and a destination t are connected, 
is known in percolation theory as the connectivity function or pair connectedness. It 
has been used in modeling epidemics or fire propagation [3 [8]. This approach is 
complementary to the effort recently devoted on complex networks, in which the network 
resilience, i.e., its robustness against link or node failures (sometimes following deliberate 
attacks) has been studied for 'scale- free' random graphs [Hj. 

Exact reliability calculations are known to be very difficult [15], except for series- 
parallel reducible graphs for which only successive simplifications {p S eries = P1P2, 
PparaiM = P1//P2 = Pi + P2 — P1P2} are needed. Even for planar graphs with identical 
edge reliabilities p and perfect nodes (i.e., p no de = 1), their algorithmic complexity has 
been classified as #P-hard [5j ITS] . Yet, the development of Internet traffic makes it 
important to assess the overall reliability of network connections, when links and nodes 
may fail. 

In this work, we show that for a network represented by an undirected graph 
G, the two-terminal reliability may be expressed as a product of transfer matrices, 
where individual edge and node reliabilities are exactly taken into account. Such a 
factorization, already observed for graph colouring polynomials [HE], 2D-percolation 
in square strips p2] or all-terminal reliability polynomials [9], UU], originates with the 
underlying algebraic structure of the graph. We apply our method to the two examples 
(K n is the complete graph with n nodes) of figure [U The .fCt-ladder describes a generic 
architecture for long-haul connections, while the i^-cylinder slightly generalizes the 
'sponge model' of width three by Seymour and Welsh [18]. When edge and node 
reliabilities are respectively equal to p and p, a unique transfer matrix is involved; 
its largest eigenvalue determines the asymptotic power-law behaviour of reliability as 
a function of the ladder length. The location of the complex zeros of Rel2(p) exhibits 
striking structure transitions as p decreases from one to zero. We illustrate the variety of 
behaviours for the above-mentioned graphs. For the sake of completeness, we finally give 
the matrix decomposition for the all-terminal reliability of the .fQ-ladder with arbitrary 
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Figure 1. Generic network architectures: (top) ^-ladder (bottom) i^-cylinder. 
Links and nodes are identified by their reliabilities: a n , b n , etc. for links, S n , T n , and 
U n for nodes. The source is So, the possible destinations are S n , T n , or U n . A missing 
link or node's reliability is simply set to zero. 

edge reliabilities (the uniform case has already been treated by Chang and Shrock [9]). 
2. Graph decomposition 

The gist of our method is to simplify the graph by removing links of the nth (last) 
elementary cell of the network, namely the edges and nodes indexed by n, a procedure 
called pivotal decomposition or deletion-contraction [5]. If the end terminal t (which 
can be regarded as perfect) is connected to node u through edge e, with respective 
reliabilities p u and p e , then 

Rel 2 (G) = (I - p e )Re\ 2 (G\e) + p ePu Re\ 2 (G ■ e) 

+ p e (l-p u )Rel 2 (GV), (1) 

where G\e and G\u are the graphs where e or u have been deleted, and G ■ e the 
graph where t and u have been merged through the 'contraction' of e; (0Q) merely sums 
probabilities of disjoint events. This procedure, along with standard series-parallel 
reductions, is repeated for the three (instead of the usual two) secondary graphs in 
order to take advantage of a structural recursivity of the graph. After a finite number 
of such reductions, we get replicas of the original graph, albeit with one less elementary 
cell and with the (n — l)th cell's edge and node reliabilities possibly renormalized by 
those of the nth cell, or set to either zero or one. In order to ensure the existence of a 
recursion relation, the graph structure must be closed under successive applications of 
(CQ); it may initially require the use of extra edges with symbolic reliabilities, so that all 
nodes of an elementary cell are connected pair-wise, even if such links do not exist in 
the graph under consideration. At this point, a recursion hypothesis is needed, giving 
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Figure 2. First step of the pivotal decomposition : the removal of edge b n . Three 
structurally identical, secondary graphs are obtained. 

for instance Rel 2 (S'o— *S n ) as a sum over specific polynomials in the reliabilities indexed 
by n; these are often obvious from the n = 2 value. Going from n — 1 to n provides 
the transfer matrix linking the prefactors of the polynomials, because Rel 2 is an affine 
function of each component reliability; the (often trivial) n = 1 case serves as the initial 
condition of the recurrence. 



3. Application to the fT 4 -ladder 



Let us first illustrate this method by calculating TZ n = Rel2(5o — > S n ) for the K4- ladder 
(top of figure [Q. Following the guidelines of the preceding section, we first consider b n 
for deletion as detailed in figure [3 Note that the three secondary graphs have essentially 
the same structure. The renewed application of leads to two families of contributions. 
The first one is a sum of 7£ n _i-like terms with prefactors, in which the 'old' a n _i, 
T n _i are renormalized by one or more of the 'new' a n , T n . The second one is a sum of 
Rekj ( Sq -^T n _i) -like terms. There is no need for coupled recursion relations for the two 
destinations S n and T n , since they are essentially identical through the permutations 
dn^e-m c n ^d n , and S n <->T n . lZ n may be expressed as the sum of five polynomials in 
a n , T n (see below). The five prefactors at step n are obtained from those at step n — 1 
by a recursion relation which translates as a 5 x 5 transfer matrix (such calculations are 
routinely performed by mathematical software). The value of 1Z\ leads to 

/ 1 \ 



Tl n = (1 0) M n M n _i ■ • ■ Mi M 







V o / 



(2) 
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where Mj's coefficients are (x = 1 — x) 



M u 


= 


(ctj + 6j ej Tj ctj 6j Tj) Si, 


(3a) 


M 12 


= 


{d; L + 6j Cj Tj <ij 6j Cj Tj) Si, 


(36) 


M 13 




dj (ij S'j + hi (xi 4" Q ej) 7j, 


(3c) 


M14 


= 


- M 44 = a , d M 42 , 


(3d) 


M 15 




- M 45 = d diM 4l , 


(3c) 


M 2 i 


= 


ifii ~\~ bi a,i Si Ci bi ai Si) Ti, 


(3i) 


M 22 




[pi 4~ bi di Si Ci bi di Si) Ti, 


(3o) 


M 23 




Ci e { Ti + bi (xi + ai di) Si Ti, 


(3h) 


M 24 


= 


- M 54 = cii ti M 52 , 


(3.0 


M 25 




- M 55 = d diM 51 , 


(3j) 


M 31 


— 


ifli bi + a j ej + &j e, 2 cij bi c^) Si Ti, 


(3*) 


M 32 


= 


- (bi Ci + bidi + Cidi-2 bi Ci di) Si T i} 


(30 


M 33 


— 


((1 - 2 bi) Xi - h {ci Ci + di di)) Si Ti, 


(3 m) 


M 34 




- M M - M 24 , 


(3n) 


M 35 


= 


- M 15 - M 25 , 


(3o) 


M 41 




Q>i bi Ci Si Ti, 


(3p) 


M 42 




bi Ci di Si Ti , 


(3g) 


M 43 




k {Xi + CiCi) SiTi, 


(3r) 


M 5 i 




&i &i Si Ti , 


(3s) 


M 52 




hi Ci di Si Ti^ 


(3t) 


M 53 




h(Xi + a i di) SiTi, 


(3u) 



with Xi = (kcldi ei + ai Ci di'cl—ai Ci di c^. In the n = case, a = 1 and c = d = e = 0. 
The five above-mentioned polynomials are actually given by the first row of M n . 
Rel 2 (So — > T n ) is given by ([2]) if the left vector is (0 1 0). We have here another 
useful instance of a product of random matrices [19J. 

The case of identical reliabilities = • ■ ■ = Ci = p (unless i = 0, see the restriction 
above) and Si = Ti = p is worth investigating, since only the nth power of a unique 
matrix needs be taken. Because of the recursion relation between successive values of TZ n , 
the generating function Q[z) = Y^=o Tin z n is a rational fraction of z. Its denominator 
T>(z) is derived from the characteristic polynomial of the transfer matrix, taken at 1/z. 
The numerator of Q(z) is then deduced from the computed first terms of the G{z) V(z) 7 s 
expansion. The final result reads: 

S(z) =\p( 1 ~PP) + ^ (4) 
Af(z) = -p{l+pp) 
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Figure 3. Location of the complex zeros of the two-terminal reliability polynomial 
T^niPt P) for n = 150 and p = l. 



\2 „5 



10p+ 13p 2 - Ap 3 -p 3 p) z 



+ (1 - p) 2 p b (2 - 4p + p 2 ) (1 -p) p 5 z\ 
V(z) =l-pp (2 + Ap p - lAp 2 p + 13 p 3 p - Ap A p 

+ 2{l-p)p 3 p 3 (2 - 7p + Ap 2 + 7p 2 p 

— I0p 3 p + 5p 4 p — p 5 p 



(5) 



4 (2-p) (l-p) 3 p 6 (1-p) p 5 z 3 . 



(6) 



Equations (J3HSD are simpler for perfect nodes, because the denominator is of degree 2 
in z\ a partial fraction decomposition provides 
1 — p 



A± 
^4 



2 

p 
2 

1 +p 
4 



5 n>0 + a+ A" + a_ AI 



2 + 4p - 14p 2 + 13p 3 - 4p 4 ± 



± 



2 + 2p+ 10p 2 - 27p 3 + 19p 4 -4p 5 

aVa 



(7) 
(8) 

(9) 



4 + 32 p 2 - 204 p 3 + 452 p 4 - 516 p 5 



+ 329p 6 -112p 7 + 16p 8 . (10) 

As n grows, 1Z n ~ a + A" : the two-terminal reliability exhibits a power-law behaviour, 
the scaling factor being A + , the eigenvalue of largest modulus. Alternatively, 1Z n ~ 
exp(— n/£), where C, = — l/ln(A + ) is the correlation length of percolation theory [7|. 

The location of the zeros of Rel2(p) in the complex plane is also worth investigating. 
The situation differs from that for chromatic [H [11] and all-terminal polynomials [9], 
because Rel2(p) is not a graph invariant. However, the node reliability p is an extra 
parameter that has a deep impact on the curves to which the zeros of Rel2(p) converge 
as n — » 00. The critical values of p at which shape transitions occur may be deduced 
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[2"U] from the three roots of T>{1/ z). The straightforward but tedious procedures used to 
determine these values, along with a few asymptotic expansions as p — > 0, are outlined 
in the Appendix (they are also applied to the i^-cylinder configuration). We limit 
ourselves to the final results in the following sections. 

A sample of the richness of behaviour is displayed in figures [3H6] for the .fQ-ladder 
and decreasing values of p. We initially observe four well-separated 'curves' that merge 
into two when p is exactly equal to 0.8 (see figure H]) and separate again. When p further 
decreases, other isolated zeros appear, as in figure El These zeros occur in pairs, the 
separation of which vanishes exponentially with n, and converge to roots of the algebraic 
equation 

= 2 + 2p + 4(3p + l)pp- (40p + 11) pp 2 

+ (45p + 4)pp 3 -20pV + 3pV. (11) 

Equation (fTTj) is obtained by ensuring that N(z) and T>{z) have a common root. The 
true limiting isolated points are such that this root is the eigenvalue of greatest complex 
modulus at the given p and p. Actually, the triplet of figure [6] appears only when 
p < p Cl pa 0.175221381869, where p Cl is a solution of (see the Appendix) 

= - 32768 - 198656 p + 3990544 p 2 - 12843528 p 3 

+ 16258037 p 4 - 6757568 p 5 - 2015436 p 6 - 575540 p 7 

+ 4636356 p 8 - 3082436 p 9 + 624640 p 10 , (12) 

whereas the associated p Cl ~ —0.604692601721 is a solution of 

= 40 - 364p + 1064 p 2 - 700 p 3 - 1946 p 4 + 4296 p 5 

- 3465 p 6 + 1074 p 7 + 146p 8 - 176 p 9 + 32p 10 . (13) 

If p Cl < p < p C2 ~ 0.406657811123 (the algebraic equation satisfied by p C2 is actually of 
degree 65 in p), only the two rightmost isolated points are present. 

The leftmost isolated point, located on the real negative axis, is asymptotically given 
by —(2 p) -1 / 3 + 25/24 + 0(p 1 / 3 ); for the other two, p must be replaced by pe ±2l7r . By 
contrast, the algebraic curves' asymptotic limit is a circle of radius (2p) _1//4 centred at 
(27/32, 0), demonstrating a different power- law behaviour with p. Finally, a third critical 
value p C3 ~ 0.49137068 also appears, for which we have not been able to find the defining 
algebraic equation satisfied by p (its degree is likely to be large); at this value, there is an 
asymptotic (anti-) crossing of the curves in the vicinity of p ~ 1.55533445 + i 0.55314582. 

4. fT 3 -cylinder 

In the second architecture of figure Sq is still the source while S n , T n , and U n are 
the three possible destinations (the last two are equivalent through a permutation of 
variables). The crucial point is to take all fi 7^ 0, because in the successive applications 
of ([1]), the merging of nodes entails a secondary graph in which S n -i an d U n -i are 
connected. As mentioned above, the dummy — with respect to the Manhattan-like strip 
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- link f n between S n and U n must therefore be present right at the start; this allows 
us to unveil the coupled recursion relations between the source and all the destinations. 
Each source-destination reliability is a sum of eight polynomials in reliabilities indexed 
by n. This could lead to 24 x 24 transfer matrices Mj. However, several rows of 
these matrices, if not identical, are linearly dependent; rearrangements of terms actually 
reduce their size to 13 x 13, even when /j = 0. 
The final result reads 

K n = v L M n M„_x • • - M x M Vfl, (14) 

where is the column vector defined by (v#) fc = <5 fcl for 1 < k < 13 (using the 
Kronecker notation: 5 k i is equal to 1 if k — I and if k ^ I), and is a row vector 
which depends on the destination: (v Sn ) k = d~i k , ( v T n ) k = o~2k, and (vu„) k = hk- The 
matrix elements are much lengthier than in (l3aH3~»j) . and are given in the Appendix for 
the sake of completeness. 

li- fi = o 

Following the procedure outlined in the preceding section, we can compute the new 
generating function. For perfect nodes, Q(So — >U n ) is given by Af/ \V1V2): 

AT = p 2 - (1 - p) p A (3 + 3p - Ap 2 ^ 

+ (l-p)V (2 + llp-3p 2 -2p 3 

+ (1 -pfp s (2 - Ap + 3p 2 + lip 3 - I3p 4 + 3p 5 

-(l-p)V° (3 + 6j9-12p 2 + 10p 3 -10p 4 + 4p 5 ) z A 



z 



Z 3 



+ {l-pfp 12 (l + 8p - p 2 - hp 3 - p 4 + p 5 ) z 5 



-(l-p) 8 p 15 [2 + hp-Ap 2 ) z 6 + (1 - p) w p 18 z 7 , (15) 
f>x = 1 - (l - p 2 ) p (l + p - p 2 ) z 

+ (l-p)V (l+p + p 2 - 2p 3 ) z 2 - (1 -p) 4 p 6 z 3 , (16) 
V 2 = 1 -p (2 + 2p + p 2 -9p 3 + 5p 4 



9 1 

+ (1 — p) p 2 [l + hp + hp 2 — Qp 3 — 15 p 4 
+ 13 p 5 + p 6 - 2p 7 ) z 2 

- (1 -p) 2 p 4 {2 + Qp + 6p 2 - 2Qp 3 + 17 p 

-18 p 5 + 27 p 6 - 16 p 7 + 3/) z 3 

+ (1 - p) V (l + 6p + 4p 2 - p 3 - 17 p 4 



A 



9p s + 3p° -2p') z q 
[2- 

l- p f p ^ z \ (17) 



[l-pfp 9 (2 + Ap + p 2 -7p 3 + 3p 4 ) z 5 
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-0.2 
-0.4 

-0.5 -0.25 0.25 0.5 0.75 1 1.25 
Re(p) 

Figure 7. Complex zeros of the two-terminal reliability polynomial for the K%- 
cylinder, /j = 0, n = 100 and p = 1. 



0.75 
0.5 
0.25 
lm(p) 
-0.25 
-0.5 
-0.75 

-0.5 0.5 1 1.5 

Re(p) 

Figure 8. Same as figure El with p = 0.6. 

When p 7^ 1, the degrees of Af, T>\, f> 2 are still 7, 3 and 6, respectively; their expressions 
are only lengthier. 

The eigenvalue of greatest modulus A max involved in the asymptotic power-law 
behaviour obeys T> 2 (lj 'A max ) = 0. The degree of the denominator leads us to expect 
that the 'width' of the network should drastically affect the size of the transfer matrices. 

The associated complex zeros are displayed for various values of p in figures [7H9J 
The overall structure is more complicated than that for the ^-ladder, but some features 
are quite similar. 

A segment of the real axis appears as a limit curve when 0.4202958 < p < 0.8092264. 
These critical values obey different criteria. Indeed, the higher one (with the associated 
critical, real p ~ 1.53039659) occurs when two complex roots of T> 2 {z) have the same 
(lowest) modulus as a real negative root of T>i(z). By contrast, the lower critical value 
0.4202958 appears when T> 2 (z) exhibits two complex roots and a real positive root with 
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2 
1 

lm(p) 
-1 
-2 



1 

Re(p) 



Figure 9. Same as figure \7\ with p = 0.1. 



the same modulus (the critical p is about 1.8363587). 

What happens when p — > ? The outermost parts of the curves tend asymptotically 

/ , \l/4 

(5 - V17) 0.684261 
to a circle of radius ^==+- — , i.e., approximately — — . The closed curve on 



fTp 



VP 



i a 



the left survives. For instance, a triple point p t goes asymptotically as ± h b, with 

VP 

a ~ 0.4610389 and b ~ —0.08457522 (a 2 is a root of a polynomial of degree 10, and 

b is a rational fraction of a). From each of these points, two curves head back to the 

ia' 

origin. One of them crosses the imaginary axis at p ~ ± — , with a' « 0.33529987 (a' 2 

\ff> 

is actually the root of a polynomial of degree 17). 



4.2. fi + 

In this case, the generating function (^(So — > U n ) is now equal to Af'/ (T>\ T>' 2 ), where 
for perfect nodes 

„2\ 



Af' 



V 



p(l + p — p 

- (2 - p) (1 - pfp 3 {l+p)(l + 3p-3p 2 )z 
+ (1 - pfp 5 (1 + 10 p + 8p 2 - 5p 3 - 2p 4 ) 2 2 

- (1 - pfp 8 (3 + Sp - 25 p 2 + 9p 3 + 4p 4 - p 5 ) 
+ (1 - p) 8 p n (1 - 2p) (3 + 3p - lp 2 + 2p 3 ) / 
-(l-p) n p 14 (l-3p + p 2 )^, 

1 - (1 - pfp (1 +p) (l+p-p 2 

+ (l-pfp i (l+p + p 2 -2p 

- (l-pfp 6 z 3 , 

1 - p (1 + 3p + Ap 2 - 23 p 3 + 23p 4 - 7p 5 ) z 



3\ z 2 



(19) 
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-0.25 0.25 0.5 0.75 1 1.25 1.5 

Re(p) 



Figure 10. Complex zeros of the two-terminal reliability polynomial for the its- 
cylinder, fi 7^ 0, n = 100 and p = 1. 




-0.5 0.5 1 1.5 2 
Re(p) 



Figure 11. Same as figure [TUl with p = 0.3. 

+ (1 - pf p 3 (1 + Qp + 2p 2 - 9p 3 

-8p 4 + 16p 5 -6p e )z 2 
- (1 - p) 4 p 6 (2 + Ap + p 2 - 15p 3 + 12 p 4 - 3p 5 ) z 3 
+ (l-p) 7 p 9 z 4 . (20) 

Note that T>' 2 is of degree 4 in z (even when p 7^ 1), so that a complete 
analytical solution for the two-terminal reliability could be obtained — but would be 
very cumbersome. 

The location of complex zeros are displayed in figures [T0l - [T2l Critical values of 
different nature occur in this case. Two isolated zeros exist as long as 0.3633122889 < 
p < 1. They do not survive in the p — > limit, in contrast with the K 4 case. For 
p w 0.3633122889, they merge with a continuous curve at p w 1.466816 ± i 0.5823927. 
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-4 -2 2 4 

Re(p) 

Figure 12. Same as figure [TU1 with p = 0.01. 



For these isolated points, the relevant p and p obey the polynomial constraint 
= - 2 (1 - 2p) (l + hp - 4p 2 ) 

- p (l + 39 p - 172p 2 + 316 p 3 - 230 p 4 + 56 p 5 ) p 
+ p 2 (12 + 42 p - 416 p 2 + 947 p 3 - 899 p 4 
+382p 5 -60p 6 ) p 2 



+ p 3 (27-45p + 33p 2 -90p 3 + 135p 4 

-77j9 5 + 15p 6 ) p 3 , (21) 



the origin of which is similar to that of (fTTl) . 

Another feature is the segment on the real axis (see figure [TT1) which occurs when 
0.016301418 < p < 0.83140245. These two critical values are actually solutions of a 
polynomial in p of degree 95, and the associated critical p's, namely 1.60638989 and 
4.56013168, are also roots of a polynomial in p of degree 95. These transitions occur 
when the equation T>' 2 {z) = has a double, real (negative) root, the opposite of which 
is also a root. 

As in the preceding subsection, the global structure expands as p — > 0. The outer 

curves tend asymptotically to a circle of radius , , ~ 0.677611 p~ 2//5 . The closed 

7 ' pi 

curve on the left also survives (see figure [12]). Here again, it crosses the imaginary axis 

• /10" 1/6 3 \ 
asymptotically at p ~ ±i I -= I . 

5. Transfer matrices for the all-terminal reliability ReU 



Nodes may be viewed as perfect in this case since the node reliabilities can be factored 
out, and simpler calculations may be done because (p]) has one less term. For the 
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.fQ- ladder, the transfer matrix is 2 x 2: 



ReU(n) = (1 0) M n ■ ■ ■ M n _ x M f J j ; (22) 

The matrix elements (Mj)^ of Mj are (a? = 1 — x) 

(Mj)n = [(aj + ei) (cj + ^) - 2 aj Q dj e*] h 

+ [(^//ei) + (cj//dj)] bi, (23) 

(Mi) 12 = aj q ei [ 1 h — H 3] &i 

(2^ Cj 

+ [(ai//ei) (ci/M)] 6i, (24) 
(Mi) 21 = [{a>i//ci) + (di//ei) - 2 (a, ej//ci di)] bi 

- (Mi)n, (25) 
(Mi) 22 = (cj + d» - 2 q dj) (^ + e» - 2 ^ <%) 6, 

- (Mi)i 2 ; (26) 

in Mo, a = 1 and c = d = e = 0. This is a special case of a multivariate Tutte 
polynomial [21] . If a« = • • • = &i = p (0 < i < n), we recover Chang and Shrock's result 
(appendix 4.2 of [9]) £/ A (z) = Af A (z)/V A (z) with 

A^^) =p + p 3 (l-p)(4-3p)z, (27) 
p A ( z ) = 1 -p 2 (l2-26p + 21p 2 -6p 3 ) 2 

+ 2p 5 (1 - p) 3 (2-p) z 2 . (28) 

The asymptotic power-law scaling factor is controlled by = (12 — 26 p + 21 p 2 — 
6p 3 + VB) with B = 144 - 640 p + 1236 p 2 - 1308 p 3 + 793 p 4 - 260 p 5 + 36 p 6 . 



6. Conclusion and perspectives 

The two-terminal reliability of undirected networks may be expressed by a product of 
transfer matrices, in which each edge and node reliability is exactly taken into account. 
This result is easily extended to the all-terminal reliability with nonuniform links, as 
well as to directed graphs. We can now go beyond series-parallel simplifications and 
look for new (wider) families of exactly solvable, meshed architectures that may be 
useful for general reliability studies (as building blocks for more complex networks), for 
the enumeration of self-avoiding walks on lattices, and for percolation with imperfect 
bonds and sites. Since the true generating function is itself a rational fraction, Pade 
approximants should provide efficient upper or lower bounds for these studies. Moreover, 
individual reliabilities can be viewed as average values of random variables. Having 
access to each edge or node allows the introduction of disorder or correlations in 
calculations. The location of complex zeros of the two-terminal reliability polynomials 
exhibits numerous structure transitions, with the possible occurrence of isolated points, 
convergence to segments of the real axis, and also an expansion from the origin as p 
goes to which obeys power-law behaviours with rational exponents which may differ 
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strongly for seemingly not too dissimilar graphs. All critical values of the node reliability 
are actually algebraic values. Finally, in a more applied perspective, let us mention that 
the failure frequency v of a given connection is another important performance index of 
networks. If equipment % with reliability pi has a failure rate \, v = J2i \v% dRe^/ dp^. 
The matrix factorization makes the calculation straightforward, since each pi appears 
in one transfer matrix only. 
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Appendix A. A few recipes on the determination of the complex zeros of 
two-variate polynomials 

Our method relies on well-known results for the zeros of recursively defined one- 
parameter polynomials [4j [9j [HJ [20]. Since we are dealing here with two-variate (p,p) 
polynomials, let us outline the procedure used to obtain the figures, the critical values, 
and the asymptotic expansions given in the text. 

Appendix A.l. Determination of the two-variate polynomial 

As shown by many published studies, the convergence of the zeros to limiting sets of 
algebraic curves is already apparent for n roughly equal to 50. To be on the safe side, 
we calculated these polynomials for n = 100 or n = 150 in order to (i) be very close 
to the asymptotic limit (ii) get a good sampling of the zeros, since — especially in the 
small-p limit — they are not uniformly distributed over the asymptotic curves when 
n — > oo (see figures [7H91 and [T2]) . 

We have calculated these polynomials using Mathematica and recursion relations 
based on the denominator of the generating function. If 

V(z) = l + b 1 (p,p)z + b 2 (p, p)z 2 + --- + b m (p, p) z m , (A.l) 

then 

r4») = Ret~ l) -b 2 ( P , p) Relt 2) b m ( P ,p) Ret' m) .(A.2) 

Knowledge of the first m polynomials deduced from the generating function allows 
the quick determination of Rel^ for a given p. p has not been kept as a parameter 
because of the explosion in the number of terms, but has been given rational values in 
order to prevent numerical errors; this gives polynomials with integral coefficients that 
may be very large (hundreds of digits sometimes). Their zeros have been obtained using 
Mathematical routine NSolve, the accuracy of which must be set accordingly (higher 
than hundreds of digits). 

Appendix A.2. Limiting curves and isolated zeros 

The zeros of recursively defined (one-parameter) polynomials mostly tend to aggregate 
close to curves such as (at least) two eigenvalues have the same modulus (the largest 
one for all the eigenvalues). Assuming that the ratio of the two eigenvalues is equal to 
e 10 , we can write 

V(z) = (1 - C e ie/2 z) (1 - C e- i9/2 z)(l-b lZ 6 m _ 2 z m ~ 2 ), (A.3) 

which must be compared with I1A.1I) . Elimination of ( and the b^s leads to a 
(polynomial) relationship between p, p and even powers of t = cos(#/2). Replacing 
t by the more practical T = cos# gives a polynomial constraint C(p, p, T) = 0. However, 
the true limiting curves are defined by only a subset of this constraint's many solutions 
for a given p and T ranging from -1 to +1, because \(\ must be the largest. In this 
context, it does no harm to investigate special points of these curves. 
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Appendix A. 2.1. Double roots of T>(z) = In our case studies the endpoints of the 
limiting curves are such that both roots are equal (9 = or equivalently T = +1): they 
are thus obtained from a subset of the solutions of C(p,p) = C(p,p,T = +1) = 0. For 
the i^-ladder with p = 1, this leads to 

= (2 - p) 2 (1 - p) 4 (1 - p + p 2 ) 2 

(4 + 32/ - 204 p 3 + 452 / - 516 / + 329/ - 112/ + 16p s ) (A.4) 

which gives the true endpoints of figure [3j -0.1175415 ± i 0.2041183, 0.7609223 ± 

i 0.5877642, 1.343654 ± i 0.3456238, and 1.512965 ± i 0.4931547 (all the solutions are 

actually roots of the polynomial of degree 8). A quicker way to find these endpoints 

dV(z) 

is to investigate when T>(z) and — - — are both equal to zero. Elimination of z from 

oz 

these two equations leads to the desired C(p,p), or more accurately, to a product of 
two-variate polynomials. Confrontation with numerical estimates of the zeros allows to 
remove spurious solutions. 

Appendix A. 2. 2. Opposite roots ofT>(z) = While they are usually not associated to 
remarkable points in the numerical plots of the complex zeros, they are nonetheless 
quite useful. Indeed, they pinpoint the limiting curves and can be obtained more 
easily because they satisfy T>(z) = and V(—z) = 0. Considering the even and odd 
components of T>{z) as functions of Z = z 2 and performing the elimination of Z gives a 
new constraint C(p,p), which is nothing but C(p,p,T = —1) = 0. This task is simpler 
because the degree of the polynomials has been divided by two (this definitely helps 
because even computer-assisted computations become ugly when the degree of T>(z) 
increases). A few real zeros may correspond to opposite roots. For the .fQ-ladder 
with p — 1, -0.2430623 and 1.527648 are indeed two such examples of intersections of 
the curves with the real axis, which may ultimately be tracked down to solutions of 
-2 - 4p+ Up 2 - 13/ + 4/ = (see figure ED . 

Appendix A. 2. 3. Real roots and segments on the real axis They are frequent features 

of the complex zeros' structure. We mentioned in the previous paragraph that algebraic 

curves may intersect the real axis at a given p, the location of which can be traced back 

to particular roots of T>{z) = 0. Whole segments of the real (positive or negative) axis 

may also occur for some graphs (see figures [8] and [TTj). It happens when, for a fixed p, two 

complex conjugate eigenvalues have the largest modulus for an extended range of real 

p's. The proper assessment of the endpoints of this segment generally requires careful, 

numerical tests of the roots of T>{£) — 0. The existence of segments of the real axis 

may be restricted to a limited range of p's or may persist down to p — > 0; it depends on 

the graph under consideration. When an algebraic curve (and its symmetrical twin with 

dC 

respect to the real axis) crosses the real axis, we have C(p, p, T) — and — Op, p, T) = 0, 

op 

because p is a double (real) root at the intersection. The elimination of T gives another 
polynomial constraint between p and p. 
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Appendix A. 2. 4- Isolated zeros, intersections with the imaginary axis, and roots of 
higher order Isolated zeros correspond to values of p and p such that the residue of 
the generating function — taken at one of the eigenvalues of largest modulus — simply 
vanishes. This implies that T>(z) and M{z) are both equal to zero. Here again, the 
elimination of z gives a constraint between p and p. In the .fQ-ladder and the K3- 
cylinder with /j 7^ 0, this leads to (fTTj) and (1211) . respectively. 

In a few cases (see figures [71 and [TO]) . algebraic curves intersect the imaginary axis, 
even for vanishing p. Noting that if p is a solution, then so is —p, we get a new constraint 
allowing the elimination of T. 

In yet other instances, sets of algebraic curves join at triple points (see figures EH9]) . 
This occurs when three roots of T>(z) = have the same modulus. These points are 
usually harder to pinpoint in practice, especially away from the real axis. 

Appendix A. 3. Critical values 

Changes — sometimes quite drastic — in the global structure of the complex zeros occur 
at particular values of p: the apparition or disappearance of real segments, isolated 
zeros, and small closed curves. These changes take place when, as p varies, different 
pairs of eigenvalues have the largest modulus. Such a situation may be described in the 
following, simplified way. Let us assume that a particular point of the complex zeros' 
structure is described by Ci(p, p) = 0. As p decreases, this feature's origin changes and 
can be traced back to another constraint (^(p, p) = 0. At the critical (p c ,p c ), both 
constraints must be satisfied. Elimination of one variable among p and p leads to the 
desired critical value. Since p is kept real, we usually eliminate p. Not surprisingly, p c 
is a root of an algebraic equation, the degree and (integral) coefficients of which may 
become quite large. For instance, let us consider the apparition of the third isolated 
zero in the .fQ-ladder configuration. Its existence is based on ( fill) , which is apparently 
satisfied for p < p c . When p is slightly smaller than p c , the isolated zero — which 
remains on the real axis — approaches the leftmost algebraic curve, which intersects 
the real axis at a point such that C(p,p,T = —1) = 0. p c and the associated p c are 
therefore defined by their obeying the following two conditions, (fTTj) and 



The elimination of either p or p leads to the defining algebraic equation for the remaining 
parameter, which can be expressed as a product of polynomials. Comparison with the 
numerical data (one can always bracket p c or p c by trial and error) allows to select the 
relevant polynomial, given in ( fl2l) and ( fl3l) . 

Obviously, the elimination procedure, which heavily relies on computer software 
(Mathematica in the present case), works best when the degrees (in the variables to 



= 4- Up + 8p 2 

+ (8p - 46 p 2 + 130p 3 - 153 p 4 + 80p 5 - 16p 6 ) p 




(A.5) 
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be eliminated) of the polynomials are not too large. A point may be worth mentioning: 
finding critical values involving only real p c and p c is usually much easier than for a real 
p c and complex conjugate p c 's, because p c is associated with T c which is seldom equal 
to ±1. We have been able to calculate the critical p c corresponding to the apparition 
of the first two isolated (complex conjugate) zeros for the i^ 4 -ladder, by considering 
the conditions C(p,p,T) = and ffTTj) . which can be decomposed in real and imaginary 
parts. This gives four equations and four parameters, namely p, T, Re(p), and lm(p). 
While it does not present any conceptual difficulty, this task may become numerically 
challenging since after each elimination procedure, the degrees in the remaining variables 
have a tendency to 'explode'. Suffice it to say that the polynomial defining this critical 
p c is of degree 65, much larger than the degree 10 exhibited by (fT2] - [T3l) . 



Appendix A. 4- Asymptotic expansions 

Our general method is to first assess numerically the expansion rate of the different 
substructures, which must behave as a negative fractional power of p (because of the 
polynomial constraints in p and p). This can be done by calculating the complex zeros 
for p equal to 10~ 3 , 10~ 6 , etc. For instance, we infer from numerical calculations that 
the isolated zeros move from the origin with an expansion rate proportional to p~ 1 ^ . 
Setting p = x P~ 1 ^ in (TTTj) gives to lowest order = 2 + 4 x 3 + 0{p 1 ^) and implies that 

X 3 = — • The leading term is therefore easily obtained, down to its prefactor (note the 
symmetry of order 3 lying at the heart of the triplet of isolated zeros). The following 
terms of the asymptotic expansion may be deduced iteratively in a straightforward way. 

As regards the sets of algebraic curves, the procedure is identical, with possibly 
different exponents. The 'best' equation to start with is obtained for opposite roots (see 

above). For instance, setting p = XP ^ 4 i n AA.5I) gives = 8% 2 — h 0{p~ 1 ^), 



implying x 4 = -■ 

Note that because the asymptotic structure is not strictly circular, the following 
terms of the expansion may depend on the argument (not only on the modulus) of the 
leading term of x- Finally, in such cases as the i^-cylinder with fi = 0, the above 
procedure gives several possible analytical solutions for x with an expansion rate in 
p -1 / 2 , with very close numerical values which makes the correct identification of the 
true prefactor quite tedious. After careful numerical tests, we finally identified the 

(S-VT7) 1/4 

expansion rate as 



V2p 

Appendix B. Transfer matrix for the f^-cy Under 

The elements rrikj of the 13 x 13 transfer matrix Mj are 
mi i = Oi Si 
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mi, 2 = Q Si Tj (6j + di /j C/j - bi d { U{) = -m 4)2 
mi,3 = e» Si C/j (/i + 6j di T t - k di Tj) 

= -m 5)3 = -77111,3 
"7,1,4 = ajTOi i2 = m 2 ,4 

mi,5 = Oj 777i )3 = m 3j n = -777 5) ii = -m llt5 

m lj6 = ai Cj Si Ti (di + 6j /, - bi di fi) = -mn j6 
mij = Q ei Si Ti C/j (6j ^ + bifi + difi-2 bi di fi) 

= "li 3) 7 = -"14,7 = -"lll,7 
777i, 8 = Oi 7771,7 = "*2,6 = "*2,8 = ™3,13 = "^8,13 

= -m 5:13 = -7777,13 = -777 9i6 = -777 9i8 = -777ii, 8 

m 21 = di Si Ti (bi + di fi Ui - b { d { fi Ui) = -7774,1 

"*2,2 = Ci Ti 

m 2 ,3 = ei Ti Ui (di + k fi Si - b t d t fi Si) 

= -m 7j3 = -7779,3 
m 2 , 5 = ^ e,i Si Ti Ui (h di + bi fi + di fi - 2 b t d t fi) 

= -m 9:5 
m 2 ,g = Ci 777 2 ,3 = m 3i9 

^772,10 = CLi (1 - bi) Ci (1 - di) Ci fi Si Ti Ui 

= m 8 ,io = m 9:12 = "7i3,io = -7774,10 
= -m 7jl0 = -7779,10 
"^3,1 = aiSiUi(fi + bidiTi-bidifiTi) 

= -"l 5) i = -"711,1 

m 3:2 = CiTiUi(di + bifiSi-bidifiSi) 

= -m 7t2 = -7779,2 
"*3,3 = ei Ui 

"^3,4 = Oi Cj Si Ti C/i (6j ^ + bi fi + di fi - 2 bi di fi) 
"13,12 = a« Ci ei Si T C/j (6* + di fi - h di fi) 

= m S ,i2 = -"T-5,12 = -"^7,12 

"^4,3 = - ei Si T C/i (6j di + bifi + difi-2 bi d t fi) 

= -"113,3 

m 4 , 4 = Oi q Si Ti (1 - 2 bi - 2 fiU + 2 bi di fi U) 
m 4 ,5 = Oi ei Si T U (di - 2b { di - b { fa - 2 d { fi 

+ 3bidi fi) = -77713,5 

"He = «i Ci ei Si T U { (~h di + f\ - 2b { fi - 2 di fi 

+ 3bidifi) = -777i3, 6 
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m i:S = di Ci Ci Si Ti Ui (di - 2 &i di + - 2b { fa 

-3difi + Abidi fi) = -mi 3>8 
m 4 , 9 = (1 - hi) Ci (1 - di) Ci fi Si Ti Ui = m 12 $ 
m 5t2 = — CiSi Ti Ui {pi di + bi fi + di fi - 2 bi di fi) 

= m llj2 

ra 5 , 4 = ai Ci Si Ti Ui (di - 2bidi - bi fi - 2di fi 

+ 3bidifi) = 77111,4 
m 5j5 = di e,i (1 - fi) Si Ui (1 - h di Ti) = m n>n 
m 5fi = aibiCi(l - di)ei(l - fi) SiTiU 

= Ci (k + di-2bi di) e, (1 - fi) Si T t Ui 
m 5 , 8 = ajm 5i 7 = m 10fi = m w ,8 = TOn,i3 
m 5t9 = - CiCi Si Ti Ui (di + kfi- h di fi) 
m 6:2 = (1 - bi) Ci (1 - di) fi Si Ti U = m 12 ,2 
™6,3 = h (1 - di) d (1 - fi) Si Ti Ui = mio.3 
m 6 , 4 = a { m 6 , 2 
m 6 , 5 = aj m 6 , 3 = m 9 ,n 
me,6 = aj c» (1 - di) e, (1 - 6, /,) 5, Ti [/, 
m 6 , 7 = Q (1 - di) d (k + fi-2bi fi) Si Ti Ui 
m 6 , 8 = a,i m 6 , 7 = m 7 , 6 = m 7 , 8 = m 9 ,i 3 
w 7j i = - diSi Ti Ui (bi di + bifi + difi-2 b { d { f { ) 

= 7719,1 = -m 8 ,i 
m 7 , 4 = di Ci Si Ti Ui (-bi di + fi - 2bi fi - 2di fi 

+ 3bidi fi) = m 9 , 4 
m 7 , 5 = di (1 - ^) ei (6i + /» — 2 6 f /i) $ T t/i 
ra 7 , 9 = Q ei 71 C/j (1 - 2 di - 2 6, fiSi + 2 bi di fi Si) 

= m 9>9 

m 7 ,u = - diCi Si Ti Ui (k di + fi- k di fi) = -m 8 ,n 
™<8,2 = -^Si Ti Ui (-k di + fi-2bifi-2di fi 

+ 3bidi fi) = mi3, 2 
^8,3 = -eiSi Ti Ui (k - 2bidi - 2bi fi - di fi 

+ 3bidi fi) = 77li3 )2 

m 8 , 4 = - diCi Si Ti U (di-2bidi + 2 fi-3 k fi 

- 4difi + 5bi difi) 
m 8 ,5 = - CLi^i Si Ti Ui (2bi + di-3bidi + fi-3 bi fi 

-2difi + Abi difi) 
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m 8 , 6 = -2aim 6J 

m 8 j = - c, e. Si Ti Ui (2 ^ + d { - 3k di + fi - 3 k 

-2difi + 4bi difi) 
m 8 ,8 = - a>iCi e, Tj C/i (-1 + 3 6* + 2 d { - 4 6j di + 3 

-hbifi-Adifi + Qbidi f t ) 
"*s,9 = - Q e» 5j Tj C/i (-d, + /» — 2 6» / f — cZ^ 

"110,1 = O* (1 - ^ (1 - /*) $ T i #i = m 12,l 
?"l0,4 = CjTOi ,l 

"iio,5 = (6j + dj - 2 6j di) (1 - /,) 5, T [/» = -ra^n 

?"l0,9 = &i Q (1 - di) Ci (1 - fi) Si Ti Ui = 77111,9 

"iio.io = Oi (1 - hi) Ci (1 - di) Ci (1 - /j) 5i ^ C/j 
wn,i2 = <»i (1 - &») Q di ei (1 - /i) 5i T C/i 
?"i2,4 = «i (1 - 6») Q (di + fi-2di fi) Si Ti Ui 
m 12 ,ii = eifn 10tl 

m 12 ,i2 = (H (1 - 6j) Q ei (1 - di /<) $ T C/i 
"^12,13 = eimi2,4 = -mi 3j i 2 
"H3,i = - a* Si T C/i (di -2bidi-bifi-2di fi 
+ 3bi difi) 

mi 3j 4 = - a, q 5i T Ui(2di - 3bidi + 2 fi - 3 &i /j 

- 5difi + Qbidi fi) 
m 13 ,g = ^ (1 - di) ei (-bi-2fi + 3b, t fi) S t Ti U t 
"^13,13 = - Ci ei 5i T C/i (-1 + 2 6j + 2 di - 3 b, L d, L + 2 /i 

-3bifi-3difi + 4bidi fi) 

All the following matrix elements are equal to zero: m^, myo, myi, mi,i 2 , "21,13, 
m 2 ,7, m 2 ,ii, m 2 ,i2, "^2,13, m 3 , 5 , m 3 , 6 , m 3j7 , ra 3 , 8 , m 3 ,i , m 4 ,n, m 4 ,i2, "Hi3, ?"5,io, ?"6,i, 

"*6,9, "*6,10, "*6,11, "^6,12, "16,13, m 7,7, "*9,7, "*10,2, "Ho,7, "*10,11, "^10,12, "^10,13, "Hl,10 
?"12,3, "112,5, "*12,6, "H2,7, "*12,8, "*12,10- 

Note that for i = 0, one must set a — 1 and c = e = 0. 
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